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The wave functions of Boson and Fermion gases are known even when the particles 
have harmonic interactions. Here we generaUse these results by solving exactly 
the N-body Schrodinger equation for potentials V that can be any function of 
the sum of the squares of the distances of the particles from one another in 3 
dimensions. For the harmonic case that function is linear in r^. Explicit N-body 
solutions are given when U{r) = —2Mh~'^V{r) = (^r~^ — C,2f~'^- Here M is the 
sum of the masses and = ^M~'^YjYjmimj{y:.i — xj)^. For general \J(r') the 
solution is given in terms of the one or two body problem with potential in 
3 dimensions. 

The degeneracies of the levels are derived for distinguishable particles, for 
Bosons of spin zero and for spin 1/2 Fermions. The latter involve significant 
combinatorial analysis which may have application to the shell model of atomic 
nuclei. 

For large the Fermionic ground state gives the binding energy of a degenerate 
white dwarf star treated as a giant atom with an N-body wave function. 

The N-body forces involved in these extraordinary N-body problems are not 
the usual sums of two body interactions, but nor are forces between quarks or 
molecules. 

Bose-Einstein condensation of particles in 3 dimensions interacting via these 
strange potentials can be treated by this method. 



1. Introduction 

Exact Boson or Fermion solutions of the quantum N-body problem in which 
every particle interacts with every other in three dimensions are very rare. They 
are almost as rare in classical mechanics although Newton solved one in Principia 
(1687) (see Cajori 1934 & Chandrasekhar 1995) and there are also some very 
special solutions such as Laplace's in which the three unequal masses describe 
ellipses about their centre of mass while at each time they make an equilateral 
triangle. However Newton's solution was for all initial conditions when the force 
on body I due to body J was of the form Fij = kmjmj{xj — x/). Newton 
reduced this problem to that of N harmonic oscillators relative to the centre of 
mass. The quantum solution is similar to the oscillator solution for solid state 
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physics. The potential energy of Newton's system is 

I <J 

Here we show that this solution may be generalised to systems in which the 
total potential energy V is any function of r. We have already explored these 
systems and their generalisations in Classical Mechanics (Lynden-Bell & Lynden- 
Bell 1999). Except for Newton's harmonic case all these systems give many-body 
forces in which the force between any two bodies is approximately linear for 
separations much less than the mean but with a coefficient that depends on that 
current mean. Despite the strange global nature of these force laws they may be 
the only non-trivial quantum many-body problems that have been solved exactly 
in 3 dimensions. Only in Newton's harmonic case do the forces reduce to simple 
pair-wise interactions. It could be argued that such global forces are unnatural, 
however, in some respects the resulting behaviour mimics that found in Nature. 
Ruth Lynden-Bell (1995, 1996) showed that these systems can be used to give 
a simple model of a phase transition that can be calculated even when is 
small. Also as we now show, such forces can be used to mimic some aspects of 
gravitation. 

For a homogeneous sphere that may pulsate in radius, a, the gravitational 
potential energy is F = — |GM^/a and the gravitational force on unit mass not 
outside a is — IvrGpx. The mean square radius of the sphere is = 3a^/5 so 

V = — (3/5)^''^ GM^/r. Now forget real gravity but adopt this form of V{r) for 
the potential of one of our extraordinary N-body problems. The force per unit 
mass on any particle at x is given by 

(o\3/2 A 
- j GMi^/r'^ = -GMx/a^ = --TrGpx 

where x stands for any one of the x/ — x. 

Thus for homogeneous spheres this choice of V{r) in our extraordinary N-body 
problem exactly mimics the effect of true gravity both for global radial pulsa- 
tions of the system and for the forces on masses within it. However if the system 
departs from homogeneity this mimicry is no longer exact. For inhomogeneous 
spherical systems the true gravitational potential energy V{r) can always be writ- 
ten —kGM'^/r with a k that depends on the radial profile. By taking that to be 
the V{r) in our extraordinary N-body problem its Virial theorem will perfectly 
mimic that of the gravitational problem but apart from the homogeneous case 
the forces on the individual particles of which the system is composed will not 
be the same in the mimic. Outside gravitational theory the concept of an effec- 
tive potential is widely used in physics and chemistry, e.g., in the shell model of 
the nucleus, in quark-quark interactions at low energy, and in modelling inter- 
molecular forces. Now we have shown that motion in these special potentials can 
be exactly calculated, they will no doubt be used as approximations in those 
applications, as well as many others. Since the N-body wave-function is known 
exactly, so is the correlation energy, but that may not be a useful general guide 
to correlation because in our systems the net force on each particle is directed 
radially toward the centre-of-mass whatever the configuration of the other parti- 
cles may be. Furthermore in many real systems the interaction between any two 
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particles is strongest when the particles are closest together while it is weakest for 
the systems discussed here. In spite of this it is possible to make systems that are 
strongly repulsive when all particles try to come close together and ones that be- 
have like gravitating systems in the sense that the overall radius obeys the Virial 
Theorem for a self-gravitating system. Even without any repulsion the exclu- 
sion principle provides support for systems of Fcrmions so with aVoc —GM'^/r 
appropriate for gravity, we find configurations of White Dwarf type. 

The N-body problems discussed here arose by direct generalisation of Newton's 
work and so skipped the developments of the intervening centuries. Wc may 
nevertheless see how they fit into those developments. Liouville (1855) showed 
that if a system of D degrees of freedom had D integrals of the motion whose 
mutual Poisson Brackets vanished, then the remaining D integrals of the motion 
could be found as quadratures. He also discovered a large class of such separable 
systems while Stackel (1890) proved his necessary and sufficient conditions for 
separability. Whittakcr (1904) gives a good description of those works while he, 
Eddington (1915) and Eisenhart (1934) helped to determine and classify such 
systems; De Zeeuw (1985) gives a good historical introduction in his thesis paper. 
Lynden-Bell (1962) and Hall (1985) developed different ideas for finding classes 
of systems with integrals or configuration invariants. Carter (1968) extended such 
results to the motion of charged particles in magnetic fields in General Relativity. 
Marshall and Wojciechowski (1988) determined those potentials in D dimensions 
for which the motion of a classical particle separates in suitable coordinates and 
the hypcr-sphcrical potential of the systems discussed here can be viewed as a 
highly degenerate member of their general Z'-dimensional ellipsoidally separable 
potentials. Evans (1990, 1991) has explored systems that are superintegrable, 
having more than D integrals for D degrees of freedom. They separate in several 
different coordinate systems and the integrals are the separation constants. This 
is the case for our hyper spherical systems (see Appendix). 

In all these works separability was achieved by changing the coordinates only. 
The idea that separability might be achievable only via canonical transforma- 
tions involving the momenta as well as the coordinates was not exploited. Thus 
Kovalcvski's top (1888) provided an unexpected new system in which the sepa- 
ration was not of the standard type. Linear soluble systems were known which 
were not of the simple separable type and Routh's (1877) Adam's prize essay on 
the Stability of Motion gives a very thorough discussion. Simple examples of both 
linear (Freeman 1966) and non-linear problems (Vandervoort 1979, Contopoulos 
&: Vandervoort 1992) that need momentum dependent transformations arose in 
stellar dynamics but it was only recently, e.g., in the work of Sklyanin (1995) that 
more general ways of looking for such systems were found. 

Meanwhile a whole body of work based around Lax (1968) operator pairs and 
the Inverse Scattering method showed that there were many previously unsus- 
pected exact solutions in both classical and quantum mechanics. This field of 
endeavour is too large to be reviewed here so the reader is referred to the re- 
view volumes Solitons (Bullough & Caudrey eds., 1980), Dynamical Systems VI 
Integrable Systems (Arnold ed., 1995) and Soliton Theory: a survey of results 
(Fordy ed., 1990). The connections between soluble models of N-body problems 
and field theory are discussed in the book edited by Bazhanov &; Burden (1995). 
In that volume quantum and classical integrable lattice models in one dimension 
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are considered by Bullough & Timonen (1995) while two dimensional models in 
statistical mechanics are discussed by Baxter (1995). 

Prominent among many exactly soluble N-body models in one dimension are 
the Toda lattice (Toda 1967, 1980) also discussed by Henon (1974) and the 
Calogero model and its generalisations. [Calogero (1971), Sutherland (1971), Bul- 
lough & Caudrey (1980), Olshanetsky k. Perelomov (1995).] In two dimensions 
Baxter's book (1982) and article (1995) contain much of interest and certain solv- 
able models in 3 dimensions have been proposed by Baxter (1986) and Bazhanov 
& Baxter (1992, 1993). Probably the most prominent soluble field theory in 2 
dimensions is that of Davcy-Stcwartson (1974), sec Anker & Freeman (1978), the 
quantum version of which is considered in Pang et al. (1990). 

Sklyanin (1995) holds out the hope that all these soluble systems may eventu- 
ally be seen as special cases of the method of Separation of Variables and produces 
some supporting evidence. 

Although the above models of interacting systems of many Fcrmions or Bosons 
can be solved exactly in one or two dimensions, the calculation that follows may be 
the only exactly solved non-trivial three dimensional N-body system yet known. 
Furthermore suitable choices of the function V{r) will allow a study of the way the 
form of interaction (albeit one of our strange global type) affects Bose-Einstein 
condensation. It should also prove possible by these methods to study the effect 
of rotation on the condensation. However this paper is solely concerned with 
the solutions of Schrodinger's equation with the correct symmetry in the wave 
function, so the statistical mechanics and Bose-Einstein condensation displayed 
by these models is not further discussed here. With the somewhat more realistic 
5 function interaction it has been studied previously for one dimensional chains 
and their continuum limits, see, e.g., Bogoliubov et al. (1994) and Bullough & 
Timonen (1998). It is of course the case that all soluble models are exceptional 
and a good example of the intricacies of non-soluble models was furnished by 
Henon (1969). 



2. N-body Solutions of Schrodinger's Equation 

Let mi be the mass of the P"^ particle and xj its position vector. Writing 
M = T,mj for the total mass, the centre of mass is given by 

X = S/X/X7 , (2.1) 

where /x/ = mj/M, which implies 

S/x/ = 1 . (2.2) 

We define 'mass weighted' coordinates relative to the centre of mass 

1/2 

rj = /ij' (x/ — x) and an associated 3N dimensional vector, r, in the space 
of all the rj by 

r = (ri,r2,r3...rAr) . 
The length of r is the mass weighted r.m.s. radius of the system since 

r2 = S/i7(x7-x)^ . (2.3) 
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This expression may be rewritten in terms of the mutual separations of the par- 
ticles since 

= ^ (x/ - x) • X/ = ^ ^ liiHj (x/ - xj) • X/ = 
I I J 

J I 

and by adding the last two expressions and halving the result 

'^^ = hYlYl - xj)^ = f^if^A^i - -^jf ■ (2-4) 

I J I < J 

In practice the r vector is constrained by the fact that the centre of mass is at 
the origin so 

S^y^rj = S^/(x/ - x) = . (2.5) 
We define three mutually orthogonal unit vectors X, Y and Z in our 3N space 

by 

X=(/xl/', 0, 0, /x^', 0,0...ti]l\ 0, O) 

Y = (O, 0, 0, ... 0, O) (2.6) 

Z=(0, 0, /i}/',0, 0, Ai2^'...0, C/x^v^') . 

Then the constraints (2.5) can be rewritten 

X-r = YT = Z-r = (2.7) 

which show that r is confined to three hyperplanes through the origin. Defining 
r = r/r then r lies on the unit 37V" sphere |rp = 1 but the X constraint confines 
it to the intersection of that sphere with the hyperplane X • r = 0, which is a 
sphere in 3iV-l space; similarly the Y constraint leaves it on the intersection of 
that SAT-l sphere with the hyperplane Y • r = 0, which is a 3A'^-2 sphere and the 
third constraint leaves it on the 3A'^-3 sphere orthogonal to X, Y and Z. 

Wc arc concerned with the N-body problems whose potential energies, V, are 
functions of the magnitude r only, so Schrodinger's equation takes the form 

The key to solving this problem lies in the right choice of coordinates. In what 
follows upper case indices run over particle labels while lower case indices run 
over coordinate-vector components. 

(a) Separation ofx 
Let Rij be an orthogonal unit 3A'" x 3N rotation matrix which rotates the 
basis vectors of our 3N space so that X, Y and Z are the last three of the new 
orthogonal basis vectors. Thus with two alternative notations and assuming the 
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summation convention over lower case indices only, 

qi = Rijrj = ^ RiJ ■ r j , (2.9) 
J 

with 

q3N-2 = X , qsN-i = Y , = Z . (2.10) 

Note 

RijRkj = Sik and ^Rij • Rfej = Sik . (2-11) 
J 

Let a run from 1 to 3N-3 (rather than from 1 to 3iV). Then the together with 
the coordinates x form a complete set of independent orthogonal coordinates for 
our system. 

We shall need the partial derivatives from (2.1) and below (2.2) 

dx/dxi = fijS , (2.12) 

drj/dy.! = {djj - l^i) S , (2.13) 
where d is the unit 3x3 matrix. The centre of mass motion will separate so our 
wave functions may be taken in the form ^ = ^(x)'0(ga . . .) so 

For Schrodinger's equation we shall need 

^ = ^ + ^^^^ ' ^'-''^ 

where the cross derivative term has vanished because ^ only involves differences 
of the coordinates xj so 'Ed'ijj/dyii = 0, i.e., ^ is independent of where the system- 
as-a- whole is. 

To evaluate the second term we need from (2.9) and (2.13) 

We check that indeed "Edtp/dxi = by summing this over / and noting that 
the two sums cancel because E/i/ = 1. We now proceed to the last term in 
Schrodinger's equation (2.15) 

Y /^/ dxi • dxi dqedqj ^ ' k i ' 

But the last term involves 

dqj [^dxij 



which is zero and 
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because R is an orthogonal matrix. Hence we have the desired expression 

a has replaced j in the final term because ip is only dependent on the first 3N-3 
of the qj, and we remember that QaQa = r'^ since X, Y and Z are all zero. 
On division by ip Schrodinger's equation now takes the form 

where q stands for the 3(iV — 1) vector qa- The equation clearly separates with the 
final three terms dependent on the qa only and the first dependent on x only, so 
it must be constant. Without loss of generality we can take the total momentum 
to be ^K. Then tp = exp(iK • x) and writing E = Et - \h'^K'^/M we find 

where q • q = r^. 

{h ) Separation of Angular Coordinates 

Equation (2.19) clearly separates again in hyper spherical polar coordinates but 
it is simplest to write them symbolically by putting q = rr and regarding r as 
independent of the angular coordinates f. We need the partial differentials 

ar/aq = r (2.20) 

dr/dq, = d/dq (q/r) = r'^iS - rr) (2.21) 
and by writing r = q/r and using (2.20) 

9/5q • r = r-^{3N - 4) . (2.22) 

We write ip = ipr{r)'ip{r) and notice that is constant on radial lines so that 
r • dip/dr = 0. Then 



and 



- ^ + (3iV - 4)r-# ^ + r-V.#i (2.24) 



dq • 3q dr"^ dr dr ■ dr 

dividing by ip Schrodinger's equation now takes the form 

[r^d'^ipr/dr'^ + {3N - Aydipr/dr] - c?r'^ + U{r)r'^ = ip'^Aiip) . (2.25) 

a is given by = —2MEfi~'^ and U{r) = —2MV{r)Ti~'^. The angular operator 
A (the hyper-angular-momentum operator) is given by 

A{'4}) = -d'^i^/dv ■ dv . (2.26) 
The angular operator A also appears in the generalised in 3N-3 dimensions 
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viz 5^/9q • 9q, so we shall first study the hyper-spherically symmetric solutions 
of d'^x/dq.- dq^ = 0. Evidently 

d^X , (3iV-4) cix_Q 
dr^ r dr 

so dx/dr = Cr~(^^~^) therefore x = Sr'"^^^'"^) where we have omitted a con- 
stant of integration which is irrelevant to our purpose. We see this is correct by 
considering the first non-trivial case in which terms other than x are involved 
which is AT = 2. Then the space of the Qa is 3 dimensional and the elementary 
solution has x = Br~^. Now our generalised knows no particular origin so 
if X is a solution for r 7^ then so is x = B\r — ro|~*^^^~^^ for r 7^ fq. We 
expand such solutions both for |r| < |ro| and for |r| > |ro| in powers and the 
coefficients of these powers are YL{r) hypersphcrical harmonics, just as they are 
in 3 dimensions. In particular the L^^^ harmonics have a power in r of either or 
j,-L-(3Af-5)_ LQQi^ijig for solutions of the form r^/(r) to our generalised = 
we see that in D = 3N — 3 dimensions 

r'^2{r'f) + iD-l)r^^{r^f)=r^Aif), 

i.e., 

L{L + D-2)f = A{f) (2.27) 

and hence the eigenvalues of A(f) are L(L + 3A^ — 5). Notice that for a two body 
problem this reduces to the L(L + 1) in 3 dimensions that we know so well. We 
shall return later to look for the degeneracies of these different eigenstates but 
for solving Schrodinger's equation the eigenvalues are sufficient. For the detailed 
separation of the 3N — 4 angular coordinates each in turn see the Appendix. 

(c) The Radial Equation 
Schrodinger's equation (2.25) now reads 

d'^^r/dr^ + (SAT - Ay-^dil^r/dr - [a^ - U{r) +L{L + 3N- 5)r-^]tPr = . (2.28) 

Now in the corresponding classical N-body problem we showed that the solution 

for the radial pulsations of the whole N-body system could be found in terms of 
the radial pulsation of the corresponding two body problem with the same U (r). 
With N = 2 we have the usual Schrodinger equation for a spherical potential 
with angular momentum i 

d^a/dr^ + 2r-^dil)2/dr - {a^ - U{r) + £{£ + l)r-^)V'2 = . 

We shall suppose that this problem has been solved and the corresponding eigen- 
values and eigenfunctions are known. Now put ?/'2 = f''^X2', then X2 obeys 

X'^ + (2/3 + 2)r-^X2 - [«' - U{r) + l{l + 1) - /3(/3 + 1)] X2 = . (2.29) 

Writing 1(1 -|- 1) - /3(/3 -M) = (£ - /3)(^ /3 1) one notices that putting 

^=\{N-2) (2.30) 

and £ = L + (3, transforms equation (2.29) into precisely (2.28). We deduce the 
surprising theorem below: 
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Theorem: If the energy levels of the usual Schrodinger equation in 3 dimen- 
sions are E{P,£), where P = 1,2,3... is the radial quantum number and i 
is the angular momentum quantum number, then the energy levels of the N- 

body problem with the 'same' 'potential' U{r) are s(P,L + |(Ar-2)). Fur- 
thermore the radial parts of their wave functions are related by ipr{P, L,r) = 



Notice that the ground states are not the same because L + |(iV — 2) can not 
be zero for N > 2. Thus the ground state of the N-body system corresponds 
to an £ of ^{N — 2) which can be a high angular momentum state of the two 
body problem. We see at once that the N-body problem will have proper energy 
levels even if strongly attractive r""^ potentials are added to U just because of 
this effective increase in the d. of the ground state. 

The theorem above is very powerful in that it enables us to use everything 
that is known about the solutions to the normal Schrodinger equation in spheri- 
cal potentials and transform it into knowledge of our N-body problems. As is well 
known, not only the generalised Kepler potential U{r) = C,r~^ — C,2^~'^ , but also 
similarly generalised square well potentials, and the U oc 5{r — ro), (5-function 
potentials as well as oscillator potentials U{r) = — ^K^r"^ — C,2f~'^ all have pretty 
solutions for the eigenvalues and eigenfunctions. All this takes over directly. How- 
ever, when C,2 is so negative that the two body problem has no ground state there 
is perhaps room for doubt as to whether we can use that potential's higher an- 
gular momentum states for our N-body problem. To allay any such doubts and 
provide a pretty illustration of the truth of our general theorem, we now solve 
the N-body equation for the generalised Kepler potential. 

{d) Generalised Kepler Problem 
We have already shown that the wave function takes the form 

■0 = exp(zk • "5tjYi{Y)il),f{r) 

where may be given a further 3A^-5 indices corresponding to various compo- 
nents of L. We must now determine ipr- 

In solving equation (2.28) we follow the standard method of solution beautifully 
laid out in the book by Pauling & Wilson (1935). Setting ar = f and C,/a = C, 
and keeping only the terms that dominate at large f we find 

d^tpr/df^ - V'r. ~ . 

So the asymptotic solutions behave as exp ± f. Only the minus sign is acceptable 
so we write ipr = ri{r) exp(— f) and derive the equation for 77 valid for all f , 

r?" + (3iV- 4)f-\' - 2?7' + I [C - (3Ar-4)] f-^-[C,2 + L{L + 3iV-5)] f-^|r? = . 

(2.31) 

We look for power series solutions of the form 

00 

p=0 
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with ao 7^ and find the recurrence relation 

+ + s + SAT - 5) - [C2 + L{L + 3iV - 5)]|ap = 

= 2{p + s) + 3(iV - 2) - C] ap-i . (2.32) 
The indicial equation has a_i = p = and yields a quadratic equation for s 

+ (3iV - 5)s - [C2 + L{L + 3Ar - 5)] = . (2.33) 

In the pure Kepler case with C2 = this yields s = L or — (L + 3 — 5) of which 
only the positive s = L solution obeys the boundary condition at the origin. For 
general ^2 the solutions are (cf (2.30)) 



s = -^(37V-5)±y[L + i(3A^-5)] + (2 = - P ± \/ + L + + C2 , 

(2.34) 

of which only that with the + sign obeys the boundary condition at r = 0. 
When < a more detailed discussion is given later. If the series for 77 does not 

terminate the asymptotic form of the recurrence relation gives Op ~ 2ap_i/p so 
7] (X e'^^ and ■i/'r is divergent at 00; so the series must terminate at ap_i say and 
in (2.32) 

C = 2(P + s) + 3(Ar-2) = 2(P + s + /3) , (2.35) 
with s given by taJiing the upper sign in (2.34) (i.e., s = L when C2 = 0). 

Remembering that ( = (/a and that = —2MEh~'^ expression (2.35) can 
be recast as the energy spectrum 

^ " ~8M (P + S + iSf " "2M T > ,0 ^1^ ■ ^^'^^^ 

^ 2P - 1 + ^(2L + 3Ar - 5)2 + 4C2 

In accordance with our theorem the energy levels with general N are given by 
putting N = 2 and then replacing L by L + /3 = L + |(A^ — 2). Of course if C2 = 
we have s = L and the theorem is then obvious from the first form. Notice that 
the theorem really applies to = —2ME/Ti^ thus we can only apply it to E 
itself if wc consider a two body problem with the same mass Af as the N-body 
problem; (^^ is also taken as unchanged since we require both problems to have 
the same U{r). However this in no way restricts us to N-body problems with M 
and C independent of N; it merely means that wc change correspondingly the 
masses M and coefficients ( in the two body problems with which we compare 
N-body problems of different N. 

Some may wish to sec the precise Schrodinger hydrogen atom spectrum with 
the correct reduced mass emerging when N = 2; to get this we must evaluate 
( in terms of Ze^. Our potential energy is F = —^Ti'^M~^C/r but this r is not 
the separation of the nucleus from the electron, R, but the mass weighted r.m.s. 
separation of them from the centre of mass. Hence r = (mnin/M'^Y^'^R where 
m and m„ are the masses of the electron and the nucleus respectively. Setting 
V = we deduce that C = 2{mmny/'^h~'^Ze^. Inserting this C into (2.36) 

along with N = 2,(2 = 0,n = P-i-L,M = m-{- m„ and putting the reduced mass 
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mmn/M = rrir the energy levels of hydrogenic atoms are given by 

rrir [Ze^f 

just as they should be. 

We now return to the question of how negative ^2 can be. Since L can be zero 

the energy of the ground state ceases to be real if ^2 < — ( ^^2~^ ) which gives 

C,2 < —\ for the familiar case N = 2. Such strongly attractive forces cause the 
particles to propagate into the nucleus and the ground state ceases to exist. It may 
be seen that the limiting case has a wave function ip r~2 near the origin which 
is easily still square intcgrable f 'il)'il)*r'^dr < oo. This is also true for the limiting 

case (,2 = — for then ip oc r~'2^^^~^^ and / ^^*r^^~^dr < oo. The limits 

are surpassed for the attractions of magnetic monopoles on the magnetic moments 
of protons and for charged monopoles attracting spinning electrons (Lynden-Bell 
& Nouri-Zonoz 1998). 



3. Level Degeneracies 

In equations (2.26) and (2.27) and the attendant discussion we showed that 
the solutions of our Schrodinger equation consisted of a hyperspherical harmonic 
in 3N — 3 dimensions times a radial function. Furthermore the hyperspherical 
harmonics of degree L in D dimensions are the coefficients of in the polynomial 
solutions of Laplace's equation in D dimensions. Thus the degeneracy of the states 
of given L and given radial quantum number P will be equal to the number of 
independent polynomial solutions of Laplace's equation of degree L inD = 3N—3 
dimensions (i.e., harmonic polynomials). To determine this number we first ask 
how many independent polynomials of degree L exist in D dimensions without 
the harmonic requirement. Each can be considered as a term of the form IljX^' 
where i runs from 1 to D and SZj = L. That number of polynomials is equal to 
the number of ways of dividing L objects into D groups where a group is allowed 
to contain no objects. If we take L units and D — 1 dividing bars then the number 
of ways of ordering them is {L + D — 1)! and if we disregard the ordering of the 
D—1 bars among themselves and the L units among themselves the number of 
ways of sorting them into groups is 

so this is the number of independent polynomials of degree L. Let fL{^a) be such 
a homogeneous polynomial of degree L in the Xa. In general 

d d 

V^/l = y] ^ -R — fh = /l-2 

where /l-2 is such a polynomial of degree L — 2 which will have G{L — 2, D) 
independent coefficients. The condition that be harmonic (V'^/l = 0) thus 
imposes G{L — 2, D) constraints on the G{L, D) free coefficients in f^. Thus the 
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number of independent harmonic polynomials of degree L in D dimensions is 
giL,D) = G{L,D)-G{L-2,D) = 

_ {L + D-3)\[{L + D-l){L + D-2)-L{L-l)] _ 

~ (i^-l)!L! " ^^-^^ 

{L + D -3)1 



{D - 2)!L! 



-(2L + D-2) 



Notice that for the familiar case D = S this gives the correct answer, 2L + 1, for 
the degeneracy of the states of given L. 

When ^2 = we have the extra degeneracy between the s,p, d, f levels of hydro- 
gen. Then a state of principal quantum number n can be obtained by combining 
a wave function of given L with a radial wave function with radial quantum num- 
ber P = n — L > 1. Thus, to find the degeneracy of states with a given n, we 
need to know the number of harmonic polynomials of degree less than or equal 
to n—1, because the n — L — l extra quanta are taken up by different radial wave 
functions. To find this number we merely add a dummy group 'o' to our sorting 
of L objects into groups and ignore the number of units, Ug, that falls into that 
group. Thus the required answer is 

gnin, D) =g{n-l,D + l) = ^i^i±^-^(2n + D - 3) . (3.3) 

For D = 3 this reduces to n^, which is the familiar degeneracy of the n^^ hy- 
drogen level before spin is considered. Thus for our problems the degeneracy of 
levels of hyper-angular-momentum L for a system of N particles is g{L, 3N — 3) 
with g given by (3.2), while for C2 = the degeneracy of the n*^ level is 
gH{n, 3N-3)= g{n - 1, 3iV - 2). 

The above degeneracies are worked out assuming that none of the particles 
are identical. In practice we are more interested in problems with N identical 
Bosons or N identical Fermions and they require wave functions with even or 
odd permutational symmetries so we now study that question. 



4. Symmetry under permutation of particles 

(a) Bosons 

For Bosons we need wave-functions that are symmetrical for the interchange of 
any two particle labels x/, xj. Both x and r have the required symmetry when the 
particles are of equal mass. Just as the magnitude of the angular momentum treats 
X, y and z symmetrically in 3 dimensions, so the magnitude of the hyperangular 
momentum L is symmetrical for particle interchange. However the components 
of the term L.y and the vector r are not symmetrical under the interchange of 
particle labels. In § 2 we found the solutions for our A/^-particle wave functions in 
the form 

i/j = exp(zk • x)YL(r)V'r(r) , (4.1) 
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where ipr depends only on the scalar magnitude L. The only term that is not 
automatically symmetrical for particle interchange is Yl^f) but even that will be 
automatically symmetrical when L = 0, because Yq is constant. Thus the ground 
state and all s-states are automatically symmetrical and are possible states for a 
system of identical Bosons. 

The states we have been discussing are not the one-particle states commonly 
considered as components of A'^-particle product states (or, for Fermions, Slater 
determinants); rather our states are themselves A^-particle states. To get a sym- 
metrical A^-particle state from one lacking that symmetry we merely add all the 
wave functions obtained by permuting the labels on the particles. But whereas 
each of our wave-functions thereby generates one boson A^-particle state, such a 
state in general comes from a number of different unsymmetrical wave-functions 
so we can no longer count the degeneracies by the arguments of § 3. However 
the arguments of § 3 connect the number of Yl functions with the number of 
polynomials that are homogeneous and both of degree L and harmonic. If we can 
count interchange symmetric polynomials independent of x which are homoge- 
neous of degree L and solutions of Laplace's equation, we have the degeneracy of 
the quantum states of hyper-angTilar-momentum L. 

Let Fl{x) = Fl{xi,X2, ■ ■ ■ xm) be a homogeneous polynomial of degree L in x 
which is symmetric under the interchange of any x/ with xj. Then 
Fl(Axi, Ax2, . . . Axjv) = A^Fi,(xi,X2,...XAr) = A'^Fl(xi, xjv, • • • X2), etc. Con- 
sider Fl{x — x) = Fl(xi — X, . . . xtv — x). It is also a homogeneous polynomial of 
degree L in x and is also symmetric, but it has the property that it is invariant 
under the transformation x/ ^ x/ + A for all / (because then x — > x-|- A). Thus 
such functions do not depend on the position of the centre of mass. However it 
can happen that Fl(x — x) is identically zero even when Fl{x) is not. For this to 
happen Fl{x) must be of the form Fl{x) = x-Fl_i(x) where Fi:_i(x) is a vector 
each of whose components is a polynomial of degree L — 1 in x which is symmet- 
ric under interchange of x/ and xj. Now let G{L) be the number of independent 
symmetric polynomials which are homogeneous of degree L in 3A^ dimensions. 
Then the number of such polynomials giving rise to non-zero Fl{x — x) will be 
G{L) less the number of free coefficients in the — x • F^-i term which we might 
expect to be 3G{L — 1). However that is not quite right because a polynomial 
with a factor xy will occur as a possibility in both the x and y components of 
Fl_i and in subtracting 3G{L — 1) we will have subtracted its number of free 
coefficients not once but twice. The same double counting will have occurred for 
polynomials with factors yz or zx so we should be subtracting not 3G{L — 1) 
but rather 3G{L — 1) — 3G{L — 2). However even that is not quite correct because 
there may be polynomials with a factor xyz. They will have been subtracted off 
three times in 3G{L — 1) but added back in three times in 3G{L — 2) so they 
will still be there and they should not be since they clearly vanish when x — x 
is written for x. Thus finally the number of independent symmetric polynomials 
which are homogeneous of degree L and independent of x is 

Gi{L) = G{L) - 3G{L - 1) + 3G(L - 2) - G(L - 3) . (4.2) 

However, we still have to impose Laplace's equation 

I 
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Now in general V^-Fi(x — x) will be a polynomial of degree L — 2 in x. However, 
since -Fl(x — x) is invariant to the transformation x ^ x + V, V^-Fl(x — x) will 
also have that property. Thus in general we may write 

V^Fl(x - x) = Fl_2(x - x) 

where Fl_2 is also symmetrical for particle label interchange since does not 

destroy that property. Thus the condition Fi_2(x — x) = will put Gi{L — 2) 
constraints on the Gi{L) free coefficients of the homogeneous L^^ degree poly- 
nomial Fl{x — x). There will be just Gi{L) — Gi{L — 2) free coefficients left in 
Fl{x — x) after imposing the harmonic condition so this is the degeneracy of the 
Yl that corresponds to the (2L + 1) with L even for the 2 boson problem. Since 
Gi is known in terms of G, we have reduced our problem to that of determining 
the number of exchange-symmetric homogeneous polynomials of degree L in 3A^ 
dimensions. This is the crux of our problem and it took us considerable thought 
to solve it. Exchange symmetry involves exchanging vectors x/ with xj, so we 
do not need symmetry in all 3A^ dimensions but only between them taken in 
triples. We shall begin our considerations with the simpler case of N bosons on 
a line with each having but one coordinate xj. We then wish to know how many 
independent exchange-symmetric polynomials there are which are homogeneous 
of degree L in dimensions. 

Let ^> = SHjx^/ be a symmetrical polynomial of degree L with A'^ factors x^/ . 

5 is the symmetrising operator which is the sum over all permutations of the 
particle labels /. Different symmetric polynomials are characterised by different 
sets of integers {£i,£2, ■ ■ ■ ^n) or partitions of the integer L into N parts, some of 
which may be zero. We construct the generating function 

CXD OO 

i?i(n,x) = ^^p(n,£)nV 

n=l t=0 

where p{N, L) is the number of partitions of L into N integers that may be zero, 

and for convenience we have defined p(0,£) = and p(n, 0) = 1. We now show 
how the theory of partitions allows us to construct the function Bi [u, x) . 

(b) Partitions of an integer L 

We learn from Abramowitz &; Stegun (1964) that 'The number of decompo- 
sitions of an integer L into integer summands without regard to order is called 
p(L)'. 

For example, five may be written 

5 = 4-M = 3-F2 = 3-M-M = 2-h2-hl = 2-M-hl-M = l-M-M-M-M, (4.3) 

so we deduce that p{5) = 7. It is easiest to work with the generating function for 
the p{L) which we call A{x). For this there is a standard result see, e.g.. Hardy 

6 Ramanujan (1918), 

OO OO 

A{x) = 1 + ^ p{L)x^ = - x^)-'^ . (4.4) 
1 e=i 
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For what follows it is essential to understand how this standard result comes 
about. To do so we rewrite the product by expansions in powers of x 

A{x) = {1 + x^ +x^ + ...) (1 + + + . . .) (l + (x^y + {x^f . . .) 
(generates units) (generates twos) (generates threes) 

(l + (xY + ...) (1 + (x^)^ + . . .) (...) ^^-^^ 

(generates fours) (generates fives) 

To see how the coefficient of in this expression is p(5) = 7 we first notice that 
we must take the 1 from all brackets after the fifth, since otherwise we would get 
too high a power of x. In the 5**^ bracket we can take the x^ term but then we 
must take the 1 from all earlier brackets. Alternatively we take the 1 in the 5**^ 
bracket. In the latter case we turn to the fourth bracket. Here we may take the 
term but that can only be coupled to the x term in the first bracket in which case 
we get the split of 5 into 4 + 1. Turning now to the third bracket and taking the 
x^ term we can take it with cither the (x^)^ bracket of the second term to yield 
3 + 2 or with the (x^)"^ term of the first bracket to yield 3 + 1 + 1. Similarly from 
the second bracket we could take the (x^)^ term with the x from the first bracket 
to give 2 + 2 + 1 or the (x^)^ term with the x^ from the first to give 2 + 1 + 1 + 1. 
Finally we could take the x^ from the first bracket to give 1 + 1 + 1 + 1 + 1. Thus 
the first bracket yields the number of ones in the sum, the second the number of 
twos, the third the number of threes, etc., and in this way the coefficient of x^ 
yields p{L) . 

However, we need the restricted partition of L into N or fewer non-zero integers 
p{N,L). These are sums of partitions pi{N,L) into exactly non-zero integers. 
Looking at (4.3) we see, for example, that pi(2, 5) = 2 = pi(3, 5). If we place a 
factor u along with each factor x^ in (4.5), then the power of u in each term will 
tell us how many parts there are in the partition generated by a particular term. 
Thus in place of A{x) we consider 

A{u, x) = {1 + UX+ {ux)'^ + {ux)^ + . . .) (1 + {ux^y + (ua;2)2 + . . .) 

(generates units) (generates twos) 



(1 + {ux^y + {ux^f + ... 

(generates threes) 



(4.6) 



Then the terms in x^ will have exactly N integers in the corresponding par- 
tition of L so (cf. 4.4) 



A{u,x) = 1 + E E^'i(^'^)^"^^ = 11(1 - 
n=ie-i e=i 



UX^) 



However, we want the number of partitions with N or fewer integers, i.e., 
p{N, L) = Y^n=i Pi(^) -^). These sums will be automatically generated if we mul- 
tiply Ahy {1 + u + + . . .) before taking the coefficient of u^x^ so 

oo oo OO' 

Si(n,x) = EEp(n,£)nV = (l+n+u^ + ...)y4(n,x) = \{{l-ux^)-^ . (4.7) 

n=l t=o e=o 
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Readers will recognise the analogy of this expression with the grand partition 
function for a gas of non-interacting Bosons. So p(N, L) can be found from the 
product as the coefficient of u^x^. Had we been intere^ed in Bosons on a line 
then p{N.^ L) would have given us the desired function G{L) but our problem is 
three dimensional. Instead of partitioning L into integers li we need it partitioned 
into integer triples {iix,iiy,'^iz) and the general term in our polynomial will 

be x^i"^ x^2'^ . . . . When we permute wc do so by exchanging 

(xi, yi, zi) as a triple with say {x2-,y2-i 2^2)- The degree of our polynomial is 



N 
1=1 



A triple (2, 1,0) corresponding to iix = 2,£iy = l,£iz = will not be permuted 
into (1,0,2) by exchanging particle labels so such triples must be regarded as 
distinct. The number of partitions of 3 into triples, such that the order within a 
triple matters but the order of the different triples does not, P3{L) with L = 3, 
is already quite a handful. Writing = for an equal total we have 




;)=(■ 



2 ^ 




8 more like that 





Ki) 


= with x,y, z components 






permuted. 




8 more like that 
= with X, y, z components = 
permuted in each triple. 




9 more like that 
(not 26 because the 
ordering of the triples 
does not count). 



Hence there are 38 partitions of 3 into triples! Now let P3{L) be the number of 
partitions of L into triples, every triple being counted as different but the different 
orderings of the same triples being regarded as the same. By analogy with (4.5) 
we consider the expression 
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A{x,y,z) = 

= (1 + .T + x2...) X (1 + (x^)! + + ...) X (1 + (x^)'^ + {x^f + ...) X ... 

X(l + y + y2...) x(l + (y2)l + (y2)2 + __,) x (1 + (j/^)! + (y3)2 + 
X (1 + ^ + ^2...) x(l + (z2)l + (z2)2 + ^^_) X (l + (z3)l + (z3)2 + ^^_) 

X (1 + (yz)^ + + ...) X (1 + (x^y)^ + (x2y)2...) 

This column x (1 + (zx)^ + (zx)^ + ...) x (1 + {y'^z^ + {y'^zf + ...) 

generates x (1 + (xy)^ + (xy)^ + ...) x (1 + (z^x) + ...) 

triples column generates x (1 + {x?z) + ...) 

M M /2\ /0\ /0\ /1\ /0\ /1\ X (1 + {y^x) + ...) 

or 1 or 2 1 1 x (1 + (z^y) + ...) 

vo; Vo; VV vo; vo; vo/ vv vv x (i + {xyz) + {xyzf + ...) 

and multiple and multiple generates triples 

combinations combinations with + £y + = 3 

thereof thereof 

By analogy with the arguments beneath (4.5) one sees how the terms of the 
third degree generate all the partitions of 3 into triples that we have just enumer- 
ated. So the coefEcient of t^ in A{t,t,t) will be p^{L), the number of partitions 
of L into triples. Furthermore A{x, y, z) may be compactly written in the form 

oo oo oo 

A{x,y,z) = n n n (i-a:w)-i . (4.8) 

p=0 q=Q r=0 

For our N boson problem we are interested not in such partitions into any 
triples but in partitions constrained to have iV or fewer triples as summands. To 
get exactly A'' summands we merely insert a « in each term as was done in (4.6). 
While to get the sum of all terms with N or less summands we have to multiply 
by (1 + n + ti^ . . .) as in (4.7). Thus putting x = y = z = t the Boson generating 
function in 3 dimensions is, writing i = p + q + r, 



oo oo oo 



^(^'*) = n n 11(1 - utp-^^^T" = n (i - ut^y'^^'^'^^'^"^ (4.9) 

p=0 9=0 r=0 e=0 

where ^{i+l){i+2) is the number of terms in the triple product withp+g+r = £. 

The required function G{L) is the coefficient of u'^t^ in B{u,t). Again if one 
only wants G{L) for L < Lmax then the infinite product can be replaced by a 
finite product up to Lmax without altering the required coefficients. 

To get Gi{L) one merely takes the coefficient of u^t^ in (1 - tfB{u,t) 
while Gi{L) -Gi{L- 2) is the coefficient of u^t^ in (1 - t^){l - tfB{u,t) = 
{l + t){l-t)'^B(u,t). 

In summary the degeneracy gsiL) of the Boson state with hyperangular 
momentum L is the coefficient of u^t^ in the expression 

(i+t)(i-t)^n(i . ("0) 

n=0 

For two particles N = 2 we find that the coefficient of u"^ is (l — 3*^) (l — t"^) ^. 
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The coefficient of in this expression is 2L + 1 when L is even and zero when L 
is odd just as it should be. This is just as in the C^^ — C^^ homonuclear diatomic 
molecule of two Bosons with every odd rotational state missing as seen in Carbon 
star spectra. 

( c ) Fermions 

The argument of § 3 relates the number of hyperspherical harmonics of degree 
L to the number of harmonic polynomials of that degree. Our experience in § 4 a 
leads us to study first the number of antisymmetric polynomials of degree L in 
dimensions. 

If Xi-X2^ . . . x^^ is a term in such a polynomial then S^/ = L. Furthermore the 
antisymmetric polynomial involving that term is the Slater determinant 



x{^ 

^1 


^^1 
X2 


^^1 
■ ■''N 


•^1 


^2 





Clearly \1 fi = ij for I ^ J then this determinant vanishes. Furthermore the 
determinant only changes sign (at most) if the 1/ arc permuted. So for any term 
that survives we may, without loss of generality, take li < £2 < £3 ■ ■ ■ < ^N- Thus 
among the Ij only li can be zero, I2 must be at least 1, £3 at least 2 and so on 
with (.^ at least — 1. For a non-zero result L = > ^{N — 1)N. 

By analogy with our study of partitions p{L) for the Boson case we now study 
partitions of L into distinct parts. Let q(L) be the number of decompositions of 
L into distinct integer summands without regard to order. Thus 5 = 4 + 1 = 3 + 2 
so that g(5) = 3. The generating functions for the q{L) is, setting q{0) = 1, 

00 00 

^g(L)a;^ = n(l + ^')- (4-11) 
e=i 

However not all decompositions of L into distinct parts lead to antisymmetric 
polynomials in N dimensions. We need a decomposition into either + (N — 1) 
unequal integers or into N unequal integers, i.e., we need the coefficient of u^x^ 
in 

00 00 00 

n (1 + ux') = E E N)u''x'^ , (4.12) 
e=o L=0 N=0 

which is the expression analogous to (4.7) of the Boson case. Again it is the 
analogue of the grand partition function for Fermions. 

The generalisation corresponding to three dimensions follows the argument for 
(4.8) and yields the generating function 

00 00 00 

Eix,y,z)= nun (1 + xW) , 

p=0q=0r=0 
p+q+r^O 

which leads analogously to (4.9) to the Fermion generating function 

F{u, t) = n n 11(1 + ^tf+^+n = n (1 + ^ . (4.14) 

p=0 q=0 r=0 e=0 
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The coefficient of u^t^ in this expression gives the number of homogeneous poly- 
nomials of degree L antisymmetric for interchanges of triples in 3A^ dimensions. 
For Fermions of spin ^ we do not need complete antisymmetry but only antisym- 
metry between particles of the same spin state. To allow for the a and f3 spin 
states being alternatives we generalise E{x, y, z) to 



oo oo oo 



n n n (1 + a^v^'^) (i + xPy^z^p) (4.15) 

p=0q=0r=0 

and in place of F{u, t) we then find F"^ 



n (l + '^^i ) ■ (4-16) 



The coefhcicnt of u^t^ in is the number of Fermionic polynomials for spin ^ 
particles with the correct antisymmetry. The arguments relating the degeneracy 
of the A?"-Fermion spin ^ wave function for the state of hyper-angular momentum 
L to this expression is the same as for the Boson case. Thus gF{L) is coefficient 
of u^t^ in the expression 

{l + t){l-tfl[{l + uty . (4.17) 

To check this formula we take the coefficient of u'^, and find for the 2 Fermion 
problem the generating function 

(l + 9i + + 3^3) (l - t'^) ~^ =Y1 + 1)*^ + 3 5^ (2^ + l)t^ , (4.18) 

icvcn iodd 

in which the coefficients of t^ will be recognised as the degeneracies of the rota- 
tional states of the hydrogen molecule with 2 protons of spin | . The even £ values 
correspond to para-hydrogcn and the odd i values to ortho-hydrogen. 

We now turn to the degeneracies in hydrogen-like potentials in which differ- 
ent L states can be degenerate. Here we need the sum of the QpiL) for all 
L < n — 1 as in equation (3.3). If we multiply our generating function for gpiL) 
hj 1 + t + t'^ + t^ + . . . = 1/(1 - t) and then pick the coefficient of u'^t"'~^ we wih 
get the required sum so the hydrogenic Fermi degeneracy is that coefficient in 

oo 

(1 + t){i - tf n(i + ui^)(^+i)(^+2) _ (4_i9) 

To find the ground state we need the first energy level n for which the coefficient 
of u^t"'~^ is non-zero. To get u'^ and no more with the lowest power of t, we 
need to use the ut^ terms rather than the 1 in all the low £ brackets since others 
come with higher powers of t. Thus if the highest £ needed is £m we require 

5]](^ + l)(^ + 2) = iV . (4.20) 

e=o 

The sum is |(£m + l) (^m+2) {£m+^) but N may not be of precisely this form for 
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integer £jn, in which case we take the lowest £m that gives 
^(•^m + l)(^m + 2)(^m + 3) > iV SO that there will be at least one term in . We 
are interested in the least power of t associated with this term. This will be 

n - 1 = + 1)(^ + 2) = lirniim + l)(^m + 2)(£^ + 3) (4.21) 

e=i 

whenever N is of the special closed shell form given by the equality. Then 

im ^ (3iV)^/3 _ 2 _ (3Ar)-i/3 ^nd n = + 1 ^ (3iV)^/3/4 as iV ^ oo. 

For large N we then find that the ground state energy behaves as 

^ = -7^ . .,..2 • (4-22) 





e 




8M 




2 




4 





To compare this energy with that of a white dwarf star we must first choose an 
appropriate value of ^ so that the potential corresponds to gravity. We showed 
in an earlier paper (Lynden-Bell &: Lynden-Bell 1999) that at high temperatures 
our systems have a Gaussian density profile at equilibrium. For such a profile the 
potential energy may be expressed in terms of the total mass Nuih and the rms 
radius at equilibrium r and is 

We therefore choose 



y=[-^) G{NmHy/r . (4.23) 



^X={i-] G{NmH?. (4.24) 



2M' V47r, 
With this choice our ground state energy level is 

where we have written M = Nrue for the mass of the degenerate electrons whose 
wave function wc have been cvahiating. This expression is of precisely the form 
found by the standard equation of state of a cold degenerate non-relativistic gas 
under its own gravity. We have, therefore, established that White Dwarfs may 
be regarded as 'superatoms' - systems with N-body wave functions which are 
solutions of Schrodinger's equation in a central potential. 



5. Conclusions 

By treating the N-body problem as a separable system in 3A^ dimensions we 
have shown that it can be solved in appropriate potentials. Marshall & Woj- 
ciechowski (1988) have given the general form of potentials that allow separability 
in many dimensions and we have concentrated on the hyper-spherical one. Of the 
many others possible, a sub-class are symmetrical for exchange of the particles. 

Whereas we have shown how to classify the wave functions by hyper-angular 
momentum, those interested in rotating systems will need to develop methods of 
classification involving the 3 dimensional angular momentum; here the methods 
of Dragt (1965) and the work by Louck k, Galbraith (1972) may prove useful. It is 
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hoped that study of that problem will throw light on Bose-Einstein condensation 
of small clusters in rotating systems. 

The fact that we are only able to treat non-relativistically degenerate white 
dwarfs suggests that an appropriate generalisation for relativistically moving par- 
ticles should be sought for systems that do not radiate gravitational waves. 

In our earlier paper on the classical mechanics of these systems we showed 
that the fundamental breathing oscillation in r separates for the far-more-general 
potentials V = Vo{r) + r^^V2(r) where V2 is an arbitrary function of all the 
coordinates which is independent of scale, A, when all the xj Xxj. While 
this is still the case in quantum mechanics the energy eigenvalues depend on 
the other motions so this separation does not by itself yield eigenvalues. How- 
ever, the possibility of solid- like and liquid-like states where Vq = \kr'^ and 

r~^V2 = k' \v I — V suggests that such systems are worthy of further 

study. 

The statistical mechanics of the systems with V = — '^^ ,r < r^, gives neg- 
ative specific heats just as those studied earlier as examples of phase transi- 
tions (Lynden-Bell & Lynden-Bell 1977). However, within one system there is no 
gravothernial catastrophe (Antonov 1962, Lynden-Bell &: Wood 1968). 

We thank Dr Dragt for sending us his papers relevant to classifying our states by symmetry and 
Drs Twarnbey and Sridhar for referring us to the Calogero model. D. Lynden-Bell is currently 
supported on a PPARC Senior Research Fellowship. 

We thank the referees for detailed comments that broadened the background of the paper 
and improved its clarity. 
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Appendix 

Separation of and commuting hyper- angular-momenta in D dimensions 
We write qi = r cos 6i, q2 = r sin 6i cos 62, Qa = r sin 9i sin 62 ■ ■ ■ sin 9a-i cos 9a, 
. . . , qo-i = r sin^i sin 6*2 • • • sin^£)_2 cos (p and qo = qo-i tan^. 
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Exact Quantum Solutions of Extraordinary N-body Problems 
The metric is given by 

ds'^ = dql + dq2 + ■ ■ ■ dqjj = 

/ D-2 \ 

= dr^ + r^ del + OidO^ + sin^ ej sin^ e^dOj + ■ ■ ■ + JJ ^^^^ ^i^^^ 
= dr^ + (^hfdel + hldOl + ... + hl_^d^'^ ^ 
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i=l 



-(D-l) 



d_ 
dr 

dr 



,D- 



-id_ 
dr 

dr 



+ r 




D-l 



k=l 




Since the hj depend only on 9i with i < j we can rewrite the angular operator A 
as follows 



A 



D-2 



sm 



_d_ 

de, 



sm 



_d_ 



,_2 9' 



(A2) 



Using (2.27), Schrodinger's equation (2.25) separates giving in D dimensions 

L{L + D-2)= Vi~^yl(^) (A3) 

multiplying (A2) or (A3) by hqj_i, -ip separates into 'il^g{6)'il^^{(f)) and the separated 
equation for gives tp^ oc exp(im(/)). If we instead multiply (A3) by /if)_2 then 
the last two terms are the only ones dependent on 9d-2 so the system again 
separates and the 0d-2 equation is Legendre's equation 



sin 



d^ 



D-2 



dOo- 



D-2 



sin^ 6*0-2 



sm9D-2 d9D-2 
with ii > \m\. 

Similarly multiplying (A2) by hj^_^ the 9d-3 equation is 



+ l)i^D-2 (A4) 



1 



d 



sin^ Od-s d9D-s 



sin^ 9d- 



dll;D-3\ iliei + l)^D-3 



de 



D-3 



-£2(^2+2)V'D-3 (A5) 



sin 9d-3 

where the eigenvalue on the right comes from (A3) applied in 4 dimensions and 
The general equation is 



sm 



d 



^-1-^ 9j d9j 



sin 



D-l-j , 



dipD- 



1-i 



sin^ 9j 



(A6) 

where the eigenvalues on the right come from (A3) applied in j + 2 dimensions. 

The operators on the left of (A4), (A5) and (A6) all have simultaneous eigen- 
values of the form + j) and commute with the energy and d/d(f). Thus we 
have D = 3N — 3 commuting operators whose eigenvalues are constants of the 
motion in D dimensions - we may of course add K and get 3Ar constants of the 
motion for our N particles. Thus our system is the quantum version of a system 
that is integrable by Liouville's theorem. 
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The above integrals of motion can all be made up of the hyper-angular-momenta 
components QaPb — QbPa- Indeed the pattern of these operators is as follows: in 
2 dimensions there is one angular momentum operator which is conserved, i.e., 
commutes with the Hamiltonian. In 3 dimensions there are two extra angular 
operators that are conserved but only the total angular momentum commutes 
with the 2D angular momentum chosen to start with. Likewise in 4 dimensions 
there are ^4x3 = 6 conserved angular momentum operators ~ 3 new ones - and 
again it is the total angular momentum that commutes with both the 3-D total 
and the 2D. In D dimensions there are ^D{D — 1) conserved angular momentum 
operators of which ^D{D — I) — ^{D — 1){D — 2) = D — 1 are new ones that did not 
occur in the D — 1 dimensional case. Of these the grand total sum of ^D{D - 1) 
squares is conserved and commutes with all former totals. Thus in D dimensions 
there are D — 1 independent mutually commuting angular momenta and these, 
together with the energy, give us the D commuting operators expected from Li- 
ouville's theorem. We note that we could have chosen any pole for our spherical 
polar coordinates; each choice gives us a different complete set of operators. 

Of course the Keplerian case V = —kM'^r~^ also separates in hyper-parabolic 
and hyper-spheroidal coordinates. For such systems we have a conserved Hamilton 
eccentricity vector (often called the Runge-Lenz vector in atomic physics) c.f. 
Lynden-Bell & Lynden-Beh (1999), equation (2.25K) 

ea = k~^M~^ ^{qaPb - QbPa)Pb - Qa/r . 
b 
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